On the origin of the pseudogap in underdoped cuprates 
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We investigate the microscopic origin of the pseudogap in the weakly doped 2D Hubbard model 
using Quantum Monte Carlo within the dynamical cluster approximation. We compare our results 
with proposed scenarios for the pseudogap. All our numerical evidence is in favor of spin-charge sep- 
aration as described in the resonating valence bond picture as the cause of the pseudogap behavior. 
Scenarios of "preformed pairs", the coupling of quasiparticles to antiferromagnetic spin-fluctuations 
and stripes are inconsistent with our results. 



Introduction The existence of a pseudogap, i.e. a 
large suppression of low-frequency spectral weight in the 
normal state of underdoped high-temperature supercon- 
ductors, is now a commonly accepted experimental fact. 
A multitude of experiments has probed the magnetic, 
thermodynamic, transport and optical properties of the 
underdoped cuprates (for an overview see Refs. Jl|, 

Some of the earliest indications of a pseudogap in the 
spin-channel were found in NMR-experiments ||, where 
the spin-susceptibility as measured by the Knight-shift 
was seen to decrease with decreasing temperature well 
above the superconducting critical temperature T c . Fur- 
ther indirect evidence of a pseudogap was found in spe- 
cific heat measurements [0. In underdoped samples, the 
electronic contribution starts to decrease with decreas- 
ing temperature in the normal state well above T c . In 
transport measurements the crossover of the linear de- 
pendence of the ab-plane resistivity in temperature to 
a stronger dependence || was taken as evidence for the 
opening of a pseudogap. However, the most direct and 
reliable measurements of the normal state pseudogap are 
angle-resolved photoemission experiments ||. A highly 
anisotropic (d x 2_ a 2-like) suppression of low-energy spec- 
tral weight persisting far above T c has been found in 
underdoped samples. This pseudogap closes at a tem- 
perature T* consistent with the crossover temperatures 
determined from other measurements. 

Theorists have responded with several scenarios ||, |j, 
|To| , pd| , |l2| , |l3| , |T4j for the pseudogap. Scenarios, based 
on the vicinity of the underdoped system to antiferromag- 
netic ordering, hold the coupling of quasiparticles to anti- 
ferromagnetic spin fluctuations responsible for the pseu- 
dogap behavior || : As a consequence of short-ranged an- 
tiferromagnetic correlations, a shadow-band forms and a 
pseudogap opens up in the density of states. 

Motivated by the fact that the pseudogap has the same 
symmetry as the superconducting gap, "preformed pair" 
scenarios associate the pseudogap with fluctuations that 
lead to (i-wave superconductivity j^j. In these scenarios 
the transition to the normal state is controlled by the 
vanishing of the supcrfluid density and not the closing of 
the superconducting gap. Due to strong pairing corre- 
lations above T c , precursors of the superconducting gap 
are seen as pseudogap in the low-energy excitations. 



Ideas involving spin-charge separation are based on the 
resonating valence bond (RVB) picture |Q. In these sce- 
narios the pseudogap is due to d-wave singlet pairing of 
spin 1/2, charge neutral fermions, called spinons. At T c 
the holons, i.e. spin charge excitations, become coher- 
ent and recombine with spinons to form electron pairs 
which renders the system superconducting. 

Further scenarios arc based on the existence of a quan- 
tum critical point (QCP) close to optimal doping that 
associate the pseudogap with a broken symmetry (see 
eg. Ref. fll]]). Other approaches that ascribe the forma- 
tion of the pseudogap to the presence of strong charge 
fluctuations are motivated by the observation of charge 
stripes in some cuprate superconductors. 

Presently, there is no consensus as to the origin of the 
pseudogap in the underdoped cuprates so numerical in- 
vestigations based on model calculations are highly de- 
sirable. Early in the history of high-T c superconductors 
it was realized that the two-dimensional (2D) Hubbard 
model in the intermediate coupling regime, i.e. where the 
Coulomb interaction is of the order of the bandwidth, or 
closely related models like the t-J-model, should capture 
the essential low-energy physics of the cuprates . 

In this report we study the pseudogap behavior of 
the 2D Hubbard model using the dynamical cluster ap- 
proximation ||, HH H, || (DCA). The DCA is a 
non-perturbative approach for the thermodynamic limit, 
which systematically incorporates the effects of nonlo- 
cal correlations to local approximations like the dynam- 
ical mean field approach (DMFA) by mapping the 
lattice onto a self-consistently embedded cluster of size 
N c . We solve the cluster problem using a combination 
of Quantum Monte Carlo (QMC) and the maximum en- 
tropy method to obtain dynamics |2C| ] . Since further in- 
sight into the pseudogap phenomenon can be obtained by 
studying the effects of lattice frustration, we also present 
results for the 2D Hubbard model with additional next 
nearest neighbor hopping t'. We show that our results 
for the pseudogap behavior in the 2D Hubbard model 
support the RVB picture and exclude other mechanisms 
for the pseudogap. 

Formalism A detailed discussion of the DCA formal- 
ism was given in previous publications JTR [l6| |l7], [H], |o| . 
The DCA is based on the assumption 12 If] that the lattice 
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self-energy is only weakly momentum dependent and can 
be approximated by a constant within each of a set of 
cells centered at a corresponding set of cluster K-points. 
The Green-functions used to calculate the self-energy are 
coarse-grained or averaged over these cells. This greatly 
reduces the complexity of the lattice problem to that of 
a periodic cluster embedded in a host which has to be 
determined self-consistently. The DCA reduces to the 
DMFA for a cluster size of one and becomes exact when 
the cluster size is equal to the system size. 

Results We present results of DCA calculations for 
the conventional 2D Hubbard model characterized by an 
overlap integral t between nearest neighbors and t' be- 
tween next nearest neighbors and an on-site Coulomb 
repulsion U. We set t = 0.25 and choose the magni- 
tude 1*7*1 < 0.5, so that the band- width W = 8t = 2, 
and study the intermediate coupling regime U = W. We 
study the initial corrections to the DMFA by setting the 
cluster size to N c = 4, the smallest cluster size which al- 
lows for d-wave pairing. We have previously shown that 
this cluster size is large enough to capture the qualita- 
tive low-energy physics of the cuprate superconductors 
Due to numerical restrictions we are currently 
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not able to perform systematic studies with increasing 
cluster size. The N c = 4 simulations presented in this 
paper should therefore be interpreted as qualitative ex- 
tended mean-field results that describe effects on short- 
length scales (within the cluster). 
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FIG. 1: Temperature-doping phase diagram of the 2D Hub- 
bard model when U = W = 2, t' = and iV c = 4. The 
critical temperature T c and Neel-temperature Tk were ob- 
tained from the divergence of the pair-field- and staggered 
spin-susceptibilities, the crossover temperature T* from the 
maximum in the uniform spin-susceptibility. 

Our results are summarized in the temperature-doping 
(T-S) phase diagram shown in Fig.0. We find regions 
of antiferromagnetism, d-wave superconductivity, pseu- 
dogap and Fermi-liquid like behaviors. In this report we 
focus on the pseudogap behavior found at low doping 8. 
The pseudogap region is characterized by a suppression 



of spin excitations p3 below a crossover temperature T* , 
defined via the maximum in the spin-susceptibility (see 
Fig.^c for t' = 0.0 and doping 8 — 0.05) when accompa- 
nied by the formation of a pseudogap in the density of 
states for temperatures T < T* (see Fig.||b). T* at zero 
doping is of the order of the magnetic exchange coupling 
J = 4t 2 /U = 0.125. To a very good approximation, T* 
is equal to the mean-field Neel temperature. These ob- 
servations indicate that short-ranged antiferromagnetic 
correlations are a prerequisite to the pseudogap. 

Thus, it is natural to explore the scenario in which the 
antiferromagnetic short-range order causes the pseudo- 
gap. To this end we study the effects of lattice frustra- 
tion, i.e. of the magnitude of the next nearest neighbor 
hopping integral t' on the pseudogap. With increasing 
|*' | the spins on the lattice become frustrated and anti- 
ferromagnetic correlations are expected to be suppressed. 
On the other hand, we expect a possible RVB spin-liquid 
state to be essentially resistant to, or even be stabilized 
by finite values of t' . In Ref. fit] it was shown that the 
ground state of the frustrated antiferromagnetic Heisen- 
berg model is almost exactly reproduced by a RVB wave 
function. In Fig. |^a we plot our results for the antiferro- 
magnetic susceptibility for different values of t'. Clearly 
antiferromagnetic correlations are suppressed as a func- 
tion of increasing magnitude of t' . 

On the other hand, the density of states plotted in 
Fig. ||b remains unchanged near the Fermi surface with 
increasing magnitude of t'. In addition, the temperature 
T* where the anomaly in the uniform magnetic suscep- 
tibility shown in Fig. 0c occurs, is essentially unaffected 
by the value of t'. Lattice frustration has little effect 
on the pseudogap. Thus we conclude that short-ranged 
antiferromagnetic correlations alone cannot be responsi- 
ble for the evolution of the pseudogap, inconsistent with 
scenarios based on the coupling of quasiparticles to an- 
tiferromagnetic spin- fluctuations Moreover, we infer 
that these results favor RVB physics over antiferromag- 
netic short-range order as the origin of the pseudogap. 

In the " preformed pairs" scenarios the pseudogap is as- 
cribed to fluctuations that lead to superconductivity at 
a lower temperature T c . In this case we would expect a 
corresponding signature in the d-wave pair-field suscep- 
tibility, i.e. an enhancement due to pairing correlations 
in the pseudogap region. In the left panel of Fig. || we 
compare the inverse d-wave pair-field susceptibility P~[ l 
of the clean system (x — 0) as a function of reduced tem- 
perature T — T c for different dopings 8 in the pseudogap 
(8 = 0.025, 0.05) and overdoped (8 = 0.20) regions. Here 
we are interested in temperatures below T* indicated by 
the arrow for 8 = 0.05. As the doping increases from 
the pseudogap to the overdoped region, P^ 1 decreases at 
fixed reduced temperature indicating that pairing corre- 
lations in the overdoped region are even more pronounced 
than in the pseudogap region. Our result in the inset 
of Figj| for the local cluster |26| equal-time d-wave pair- 
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FIG. 2: (a) The antiferromagnetic spin-susceptibility versus 
temperature, (b) the density of states at fixed temperature 
T = 0.022 in the pseudogap regime and (c) the uniform mag- 
netic susceptibility for S = 0.05 and for different values of 
the next-nearest neighbor hopping amplitude t' . Frustration 
suppresses antiferromagnetic correlations but does not affect 
the pseudogap. 



field susceptibility P l d ocal = ^ J2 ijk Sij to^c^cj.^), 
where gij is the Fourier-transform of cos k x — cos k y , 
supports this finding. This quantity does not contain 
phase fluctuations and thus can be taken as direct evi- 
dence against strong pairing correlations at underdoping. 
These results lead us to eliminate the "preformed-pairs" 
scenarios as an origin for the pseudogap in our results. 

Among other scenarios for the pseudogap include the 
idea of QCP, stripes, and spin-charge separation. In most 
scenarios, the QCP exists at the zero temperature termi- 
nus of T* versus doping. However, near this point, we 
find no evidence for a QCP in our simulations; i.e. none of 
the many susceptibilities that we measure are enhanced 
as we approach this point. We do not see any signa- 
tures of stripes in the charge susceptibility (not shown) . 
However, stripes might be important for larger clusters. 

As already discussed, the results for the frustrated sys- 
tem are consistent with a RVB, i.e. spin-charge separated 
state. In order to address the question of spin-charge sep- 
aration unambiguously, we would have to calculate the 
lattice dynamic spin- and charge- structure factors for a 
range of momenta and extract the spin and charge dis- 
persions. However, these calculations, although formally 
possible, are presently numerically too expensive to per- 
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FIG. 3: Left: The inverse pairfield-susceptibility (inset: local 
cluster equal-time d-wave pairfield-susceptibility) versus re- 
duced temperature for various dopings for t' = 0. Right: The 
density of states for the same doping levels at T = 0.022. 



form. Nevertheless, we believe that we can gain insight 
into this question by studying the behavior of the corre- 
sponding quantities calculated on the cluster only[p6[. 
Preliminary results (not shown here) for the dynamic 
cluster spin- and charge susceptibilities at low doping 
are consistent with spin-charge separation: The spin- 
susceptibility, dominated by fluctuations at Q = {u,ir), 
becomes suppressed at low frequencies with decreasing 
temperature well above T c . In the charge-susceptibility, 
which is strongest at q = (0, 0), weight builds at low fre- 
quencies as the temperature is lowered and no sign of a 
pseudogap is seen. 

The results shown in Fig. [| for the static uniform lat- 
tice spin- (circles) and charge-susceptibilities (squares) 
versus temperature at underdoping (8 = 0.05, left) and 
overdoping (5 = 0.30, right) support this picture. As 
already discussed, the spin-susceptibility Xs in the un- 
derdoped system shows an anomaly at a temperature 
T* below which spin-excitations are suppressed. The 
charge-susceptibility \ c in the underdoped system how- 
ever displays qualitatively different behavior: When the 
spin degrees of freedom become suppressed at T % T* , \c 
starts to rise with decreasing temperature. This clearly 
indicates that coherence in the charge channel starts 
to set in with decreasing temperature whereas the spin 
excitations become suppressed. These results are con- 
sistent with recent angle-resolved photoemission exper- 
iments [ p7| where superconducting order in the under- 
doped system was found to be accompanied by an emerg- 
ing quasiparticle coherence as seen in the spectra. 

Thus, our results for the underdoped region can be in- 
terpreted within an RVB, spin-charge separated picture, 
where the development of the pseudogap originates in the 
pairing of the spin degrees of freedom and superconduc- 
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tivity is driven by coherence of the charge excitations. 
In the overdoped region (right panel in Fig. ||) our re- 
sults show that the system becomes more conventional, 
i.e. spin and charge excitations behave in a similar way. 
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FIG. 4: Temperature dependence of the uniform spin and 
charge susceptibility in the underdoped (left) and overdoped 
(right) regime when U = W = 2, t' = and N c = 4. 

Summary We have investigated the microscopic ori- 
gin of the pscudogap in DCA simulations of the under- 
doped 2D Hubbard model. We find that short-ranged 
correlations alone, although necessary for the pseudogap 
to emerge, is not sufficient to describe its origins. Simi- 
larly, although stripe fluctuations may emerge for larger 
clusters, no indications were found in our calculations, 
and therefore stripes are not the origin of the pseudogap. 
Likewise, we find no indication for strong pairing fluctu- 
ations for temperatures T c < T < T* , in disagreement 
with pre-formed pairing scenarios. Finally, although we 
cannot exclude a quantum critical point at finite doping 
as an origin, we find no evidence for it. Instead, all of our 
numerical evidence points to spin-charge separation as 
described by Anderson's RVB theory as the microscopic 
origin of the pseudogap. In particular, in the pseudo- 
gap region, spin degrees of freedom are suppressed while 
charge excitations become more coherent with decreas- 
ing temperature. At overdoping where no pseudogap is 
found, spin and charge behave qualitatively similar. 
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